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C^ Abstract. We study a recently proposed scenario for the early universe: Subluminal Galilean Gen- 

^=^ esis. We prove that without any other matter present in the spatially flat Friedmann universe, the 

-^ perturbations of the Galileon scalar field propagate with a speed at most equal to the speed of light. 

This proof applies to all cosmological solutions — to the whole phase space. However, in a more realistic 
situation, when one includes any matter which is not directly coupled to the Galileon, there always 
^^ exists a region of phase space where these perturbations propagate superluminally, indeed with ar- 

^-H hitrarily high speed. We illustrate our analytic proof with numerical computations. We discuss the 

implications of this result for the possible UV completion of the model. 
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1 Introduction 

The theory of generation of quantum cosmological perturbations [1-10] during early-universe inflation 
[11-15] received additional strong support from the recent results of the Planck collaboration [16]. 
Inflationary spacetimes usually feature a singularity or in other words strong quantum gravity stage 
far enough back in their past [17]. For a universe with no spatial curvature, the only known method 
of avoiding a strong quantum gravity stage is to posit a period in the evolution of the universe 
where the Null Energy Condition (NEC)^ is violated. However, until recently, such classical and 
local constructions in the context of standard general relativity were always plagued by some sort of 
pathological instabilities — either ghosts or gradient instabilities or both, see e.g. [18].^ The situation 
has changed with the rediscovery [23-30] of more general single scalar-field theories [31], which contain 
higher derivatives of the scalar field in the action but still have only one scalar degree of freedom. 
These theories were explicitly demonstrated to be able to violate the NEC with neither ghost nor 
gradient instabilities in at least a part of their phase space, "^ see e.g. [28, 33-35]. Moreover, it was 
shown in [36] that these theories can realise a bounce and avoid a strong quantum-gravity regime but 
only at the price of a strongly coupled regime of these non-renormalizable theories in the initial, now 
collapsing, universe. 

An interesting and unusual scenario for the initial stages of the universe was proposed in [34] : 
Galilean Genesis. There, the universe starts from the Minkowski spacetime which corresponds to an 
unstable and singular configuration of the Galilean scalar field with zero energy density. Small depar- 
tures from this original state violate the NEC. Thus starting from this configuration the universe can 
expand by generating the energy density during its evolution. In fact this evolution is just a super- 
inflationary stage with a very small acceleration.'* Nearly scale-invariant cosmological perturbations 
are supposed to be generated by introducing an extra spectator scalar field coupled to the Galilean in 
such a way that effective metric reproduces quasi de Sitter. Eventually the Galilean should disappear 
by transmitting its energy into radiation, before the universe runs into a Big Rip singularity [39, 40]. 
This may happen through a reheating mechanism proposed in [41]. However, it is not clear whether 
one can really get rid of the Galilean, so that it does not spoil the late evolution of the universe, 
because some of the Galilean's configurations have negative energy densities while on some other the 
NEC is broken so that the energy density in them will still grow with time even after the end of 



-"^This condition states that T'^'^rinrti, > 0, where T^"^ is an energy-momentum tensor n'' represents all null vectors. 

^One may try to avoid these problems by adding higher derivatives directly to the EFT for perturbations see e.g. 
[19-22], This approach may be useful for a better systematisation of the perturbative expansion but cannot elucidate 
the behaviour of the cosmological background. 

^However, it does not mean that there are no instabilities at all. Indeed, there are good physical arguments [32] that 
there are more subtle and less dangerous run-away instabilities which are still present in these theories. 

^This is also the case for other Genesis scenarios [37, 38]. 



reheating. In [34], only the particular Genesis solution was studied, while bouncing trajectories in 
the model were first investigated in [42]. The full cosmological phase space was investigated in [36], 
demonstrating that, without a reheating mechanism, all trajectories eventually evolve to the configur- 
ations where the sound speed becomes zero, so that the original effective field theory (EFT) becomes 
infinitely strongly coupled. This strong-coupling regime occurs before the would-be Big Rip future 
singularity. This conclusion also applies to the Genesis solution which turns out to be a separatrix in 
the cosmological phase space of the Galilean Genesis. 

Perturbations propagating around the Genesis solution initially do so with a sound speed equal 
to that of light, slowing down as the expansion becomes significant, and eventually vanishing. In 
these higher-derivative scalar theories, superluminal propagation is common, see e.g. [43]. While its 
presence does not necessarily signify the existence of causal paradoxes or problems on the level of 
EFT see e.g. [44-48], it does mean that any EFT where a connected part of the phase space features 
superluminality cannot have an ultra-violet (UV) completion which is local, Lorentz invariant and 
corresponds to some weakly coupled heavier particles [43] . Thus such theories cannot have a standard 
Wilsonian UV completion, but may still be completed in a non-standard way, via e.g. classicalization 
introduced in [49, 50] for such derivatively coupled scalars. Moreover, it was argued [51-53] that 
classicalization may only work provided that superluminality is possible. 

To avoid the superluminality of perturbations and give a hope for a Wilsonian UV completion, 
in [54], Suhluminal Galilean Genesis was introduced. An additional parameter in the action ensures 
that the sound speed on the Genesis trajectory can be arbitrarily small, while the substantial features 
of the original Galilean Genesis are preserved. This should imply that in other parts of the phase 
space, at least those close to the Genesis trajectory, the sound speed would be similarly reduced, in 
principle removing any such potential superluminality. 

In this paper, we extend our phase-space analysis presented in [36] to the Suhluminal Galilean 
Genesis model. 

First, in section 4 we analytically prove that any configuration of a spatially- fiat Friedmann 
universe filled solely by the Galilean with the action and parameters from Suhluminal Galilean Genesis 
[54] has a subluminal sound speed. In fact, this is true even of the original Galilean Genesis model 
[34]. 

However, as a result of the imperfect nature of the fiuid [55] described by the scalar field or 
kinetic mixing / braiding of the scalar with the metric, the local properties of the Galilean do not just 
depend on the local scalar-field variables, but also on the local Ricci tensor and, through the Einstein 
equations, on all other external matter species present locally. In particular, this implies that the 
sound speed of the Galilean changes in the presence of other matter species. Using this unusual 
property, we show analytically in section 5 that for all admissible parameter values of the Suhluminal 
Galilean Genesis model, there are always cosmological configurations which feature superluminal 
propagation given the addition of an appropriate external positive energy density with a normal 
equation of state Q < w < \. Moreover, perturbations around these configurations are not ghosts and 
some of these configurations correspond to the same values of the local scalar-field variables available 
in the phase space of the Suhluminal Galilean Genesis without any external matter. In particular, 
these configurations can be obtained by a continuous transformation from the Genesis trajectory. We 
illustrate our analytic proofs with numerical computations. 

Thus as such, none of the Suhluminal Galilean Genesis models is free of superluminality in the 
presence of external matter and by the arguments of [43] cannot enjoy a standard Wilsonian UV 
completion. 



2 Model, Results Summary and Discussion 

We study the class of models for the Galileon scalar-field tt, minimally coupled to gravity g^^i,, the 
dynamics of which are given by the action introduced in [34, 54] ^ 



-fe^'' {dirf + 7 ^ {dn)^ + 7 {dirf On 



(2.1) 



where D = ,9'"'V^V,y, with V^ representing the covariant derivative while 

P^l + a, 
is a dimcnsionlcss parameter bound to be in the range 1 < /3 < 4, and another dimcnsionless parameter 

7 = ({) » 1 , (2.2) 

is a combination of constants / and A of dimension one. The positive parameter a was introduced in 
[54] and represents the only difference between the system (2.1) and the so-called Conformal Galileon 
used to drive the Galilean Genesis scenario in [34], where a = 0, and which was also studied in context 
of a bouncing cosmology in e.g. [36, 42, 56]. The reason for the introduction of the positive a in [54] 
is the "robust subluminality" around the Genesis solution. We will mostly use /3 to simplify formulae. 
It is useful to perform the following field redefinition 



vr - In ( ^ ) , (2.3) 



so that one obtains a form of the action where the scalar field (j) has standard dimensions, 

\4 /o^\2 



5"^ = / d'^x^/^t 



2 , _ I P-2\ {dcj^r , {dcby n</) 



W)^ + 7(^)^^+7 



(2.4) 



and the action has just two dimcnsionless parameters /3 and 7. In particular, from this form of the 
action, it is clear that no single physical observable related purely to the evolution of this system can 
depend on A and / separately but only on their ratio given by 7. However, if there are non-minimal 
but non-universal direct couplings to some external fields / matter through 

9T'" = ^^^9^..= {J^ 9^., (2.5) 

as was proposed in [34] for a spectator field a, then / plays a role of the coupling constant to this type 
of external matter, see [25]. This coupling is not only needed to generate fiuctuations as in [34] but 
also to reheat and possibly exit from the super-inflationary regime [41] . Note that the coupling cannot 
be universal for all matter, since the model (2.1) would then just correspond to inflation written in 
a conformally rescaled metric. Therefore there should be at least some matter species which do not 
couple to the Galileon directly through (2.5) but only through gravity. Further we note that the 
system simplifies for 13 — 2. 

The dimcnsionless constant 7 can be moved in front of the whole action if one rescales the 
coordinates as 

x^ -^ 7^/^ y'' , so that 50 ^ 7 50 . (2.6) 

Note that the Einstein-Hilbert action scales in the same way S'eh — > 7'5'eh. Hence the dynamics of 
the Galileon field 4> interacting with gravity is completely independent of 7 in the absence of other 

^Contrary to [34, 54] we use the signature convention (+,—,—,—), so that the sign in front of the first term is 
"wrong" — it corresponds to a ghost. 



fields. In a realistic universe there should be other matter species, e.g. radiation or cold dark matter. 
If we assume that those are not coupled to the Galileon directly but only through gravity one should 
rescale their energy density (and pressure) by 

P^P/I, (2.7) 

so that, after this rescaling, 7 disappears from all equations. Therefore the rescaled energy density is 
allowed to be larger than unity but should be parametrically smaller than 7 in order to avoid physical 
transplanckian energy densities. We define two rescaled variables, which we will use to describe the 
system, 

TO = TT 7I/2 , and h = H 7I/2 . (2.8) 

Then, following the method introduced in our paper [36], it is illustrative to describe dynamics in 
phase space {ra,h) rather than (tt,™) by solving the Friedmann equation with respect to e^^: 



„27r 



e 



3to4/3 + I2hm^ - Qh 
2^ 



2 



>0. (2.9) 



Only those regions in (m, ft.)-phase space where the r.h.s. of (2.9) is positive are dynamically accessible 
to the system. 

This method allows us to present the first result immediately: we have plotted the phase space in 
the (rn, h) coordinates in Fig. 1. The red regions are dynamically inaccessible, but there are no ghosty 
regions in the phase space at all. The Galilean Genesis trajectory is marked as the thick salmon line 
originating at the origin. The regions of superluminality are coloured in deep blue. We can clearly see 
that for /3 > 1 they never occur in the phase space. Indeed there is no superluminality anywhere in 
the phase space even in the original Galilean Genesis scenario of [34]. We have proven these results 
analytically in section 4. 

The situation changes upon adding external matter with positive energy density p, coupled to 
the scalar tt only through gravity. In Fig. 2, we have plotted a two-dimensional slice through the now 
three-dimensional phase space {m,h,p), keeping p fixed. Adding external energy density opens up a 
ghosty region surrounding the origin of the (m, h) axes, which we have plotted in yellow. The boundary 
of this region is a pressure singularity at which the sound speed diverges. This neighbourhood of this 
boundary can emit trajectories and can be approached by trajectories. Therefore there is always a 
superluminal region, containing arbitrarily high sound speeds, surrounding this ghosty region. This 
is true for arbitrary small positive external energy densities p and for all values of /3 relevant for the 
Subluminal Galilean Genesis scenario. Adding external energy density allows for new cosmological 
configurations of the scalar field (7r,7r), opening up the phase space in this sense. The addition 
of positive energy density with a normal equation of state can sometimes transform a cosmological 
configuration of the scalar field, into a configuration without ghosts but with superluminal propagation 
of perturbations. We prove this analytically in section 5. 

Further, we would like to mention that, similarly to the Ghosts Gondensate [57] or the simplest k- 
inflation [58], the systems described by (2.1) or (2.4) do not possess any physically acceptable Lorentz- 
invariant vacuum. This breakdown is acceptable for a theory effectively describing a hydrodynamics of 
an (im)perfect fiuid. However, as a field theory it is far beyond the Standard Model — the QFT which 
describes the current state of knowledge in particle physics. The appearance of ghosts around X = 
separates the phase space into two disconnected sets which correspond to two different EFTs — one for 
configurations with spacelike derivatives and another one for the time like case. Note that a smooth 
dynamical evolution from one set to another, i.e. from timelike to spacelike derivatives or vice versa, 
is physically impossible, because of the pressure singularity preventing these transitions [28, 36, 55] 
in these theories without a healthy Lorentz-invariant vacuum. This is obviously the case for the 
standard representation of fiuid dynamics through fields, see e.g. [59]. Close to normal vacuum, where 
the energy density of the fiuid and particle density are zero, the fluid description should break down 
in any case and one is forced to use kinetic theory. In field-theoretical language this breakdown would 
correspond to some pathology of the EFT — to an infinitely large strong-coupling scale or ghosts or 
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Figure 1. Phase portrait for the Subluminal Galilean Genesis system (4.4) for various values of the 
parameter /? = 1 + a with no external matter present. The axes are rescaled coordinates h = {f/A)^'^H, 
m = {f/A)^'^-k. The deep red region is dynamically inaccessible, the blue corresponds to configurations 
where the sound speed is superluminal, the sandy region — to those configurations where c^ < and gradient 
instabilities are present. Healthy trajectories evolve through the white regions and have been illustrated with 
flow lines. The Galilean Genesis trajectory discussed in [34, 54] is a separatrix and has been explicitly singled 
out and marked as the thick salmon-coloured line. Nowhere in this phase space is there a region where the 
perturbations are ghosty. 

For /3 > i, the accessible phase space never has any superluminality anywhere. This includes the ori- 
ginal Galilean Genesis model with ,5 = f. As ,5 increases, the sound speeds are reduced and eventually for 
/3 = 4 the sandy region with imaginary sound speed reaches to the origin of the phase space which should 
correspond to a singular Minkowski space. For /3 < 1, an accessible superluminal region appears at the origin 
and grows in extent with decreasing /3. 




Figure 2. A slice through the phase space of the Subluminal Galilean Genesis model with /3 = 2 in 
the presence of external radiation (w = 1/3) with rescaled energy density p = 0.3. The axes are rescaled 
coordinates h = {f/A)^'^H, m = (//A)^"7r. The whole phase space is {m,h,p). The deep red regions are 
dynamically inaccessible; inside the yellow region the perturbations of tt are ghosts. The sandy region has 
Cg < while in the blue the sound speed is superluminal. The sound speed diverges as the boundary of the 
yellow ghosty region. 



The blue superluminal region to the top right of the yellow ghosty region is accessible to trajectories 
and has arbitrarily high sound speed close to the ghosty region. It appears for arbitrarily small values of p 
and even at ft = 0. We demonstrate this analytically in section 5. 



unacceptable gradient instability. This property is crucial for the EFTs which allow for superluminal 
propagation. Indeed, if such EFT can only exist for timclikc gradients of some scalar field, then one 
can always declare this field to be the time coordinate so that the resulting spacetime is stably causal 
and any causal paradoxes are impossible by definition, see e.g. [46]. 

On the other hand, exactly the absence of a physically acceptable Lorentz-invariant vacuum in 
the context of an EFT implies that one cannot use this EFT to describe a transition from inflation 
(or from some Genesis stage of the early universe as in [34, 54]) to standard hot Big Bang cosmology. 

Properties of a possible UV completion of an EFT should depend on all states / configurations 
this EFT is able to describe. It is not enough to look at a given state or given trajectory and its 
small neighbourhood. Indeed, why should the EFT be inapplicable further away if nothing prevents 
it from operating there? Only in the case when the phase space (or Hilbert space) has different 
unconnected regions can one claim that there are different EFT's describing these regions separately. 
These regions can be separated by states or configurations with ghosts or gradient instabilities or 
some classical singularities or infinitely strong coupling. For example, this is the case for perfect 
fluids where the same Lagrangian could be formally used for spacelike gradients. In some of these 
configurations with spacelike gradients, superluminal propagation may be possible. This, however, 
would not at all imply that one cannot find weakly coupled particles building the fiuid corresponding 
to the timelike gradients for which there is no superluminal propagation. 

In this paper we have only studied external matter which is not coupled to the Galilean. Indeed, 
this matter should unavoidably be present in the model. However, we do not think that a non-minimal 
direct coupling (2.5) would change the situation. 

Another potential way to avoid our conclusion is to prove that the energy density p^ (</), m) 
needed to induce superluminality around a configuration ((/), m) brings the EFT away from its region of 
validity. However, this is hard to check, since the strong coupling scale fi (for wave vectors) generically 
depends on all cosmological phase space coordinates, so that fx = fj,{(j), m, p). On dimensional grounds, 
one would expect that the limiting external energy density for a given ((/), m) should be provided by the 
solution of the equation p = Cs {4>, m, p) p"^ (0, m, p), where Cg is the sound speed. Unfortunately there 
is no self-consistent derivation of /^ as a function of an arbitrary cosmological configuration (0, m^ p). 
This analysis is definitely interesting and important but involves the calculation of the cubic action for 
the cosmological perturbations in this general non-slow-roll setup. This task goes beyond the scope 
of our paper. Moreover, it follows from our analysis that the corresponding ps {(j), m) is not necessary 
parametrically large. If the addition of such an external p can invalidate the EFT it could imply that 
reheating might also lie outside of the region of validity of this EFT. 

It would also be interesting to investigate whether one could modify the theory in such a way 
that superluminal propagation would be impossible on all configurations with timelike gradients in 
the presence of any matter with a normal equation of state. In such a case, the superluminality could 
only occur for some unrealistic external matter. However, given the properties above, this could be 
impossible and would definitely be a cumbersome and ambitious task. But only such strong property 
could give a hope for a standard Wilsonian UV completion without a need for something less standard 
like classicalization [49, 50]. 

We expect that the situation with UV completion is typical for all recently rediscovered scalar field 
theories able to violate the NEC. 

3 Some properties of models with kinetic gravity braiding 

Both the versions of the action for Galilean Genesis, (2.1) and (2.4), belong to the same class ^ 

S^= f d^x^^ [K {if, X) + G ((p, X) n^] , (3.1) 

®This class of models was introduced in [28] as kinetic gravity braiding and then slightly later in [35]. 



where K and G are arbitrary functions of ip and where 



X^-g^^•'^^<p^,^. (3.2) 



We are going to work with both actions (2.1) and (2.4), hence, to make the paper self-contained, we 
will list the main dynamical equations and stability criteria valid for a general theory of type (3.1). 
After that one can use these equations by substituting ip ^ n or (f) and the corresponding K and G. 
All these equations can be found with derivations in this form in [28, 55]. In a cosmological setup, for 
these systems, it is convenient to use hydrodynamical notation and analogy introduced in [55]. 
The diffusivity measuring the kinetic mixing / braiding of ip with the metric is 

K = 2XG^x , (3.3) 

where {) -^ = d{) /dX. The effective mass / chemical potential is 

m = dt'p = p ■ (3.4) 

The pressure of this scalar imperfect fluid is 

r ^K-m^G,^- Km, (3.5) 

while the density of shift charges (which are not conserved for the systems studied in this paper) is 

n = K^m — 2mG^^ + 3Hk , (3.6) 

where H is the Hubble parameter. The corresponding energy density is 

£ = ran — V — nm = e + 3HmK , (3.7) 

where 

e = m{K,^-mG^^)-K, (3.8) 

is a "perfect part" of the energy density which is independent of the expansion or Hubble parameter. 
In the spatially flat Friedmann universe we have ^ 

H^^^i£ + p)^KmH+^{s + p) , (3.9) 

where p is the energy density in external matter species. We have assumed that these species are 
not coupled to the Galileon directly, so that the corresponding equation of motion is the continuity 
equation p + 3H (p + p) = 0, where p is the their pressure. The second Friedmann equation is 

H = --{£ + p + r+p) = -{Km-nm-{p + p)) . (3.10) 



The equation of motion for the scalar p is 



1 

2'-- / ■ -,v 2' 



m D + 3n ( iJ — -nm ) + £^^ = -k{p + p) , (3.11) 



where 

D^^-^ + iHn^ + -K^ . (3.12) 

TO ' 2 

The partial derivatives are all taken at constant Hubble parameter H . The perturbations are not 
ghosts provided D > Q. Finally the sound speed is given by the formula 

2 _ n + k+ k{H - Km/2) 

Dm • ^■^■^■^> 



'^Throughout the paper we use the reduced Planck units where Mpi = (SttGn) ' = 1. 



Fluids with Kinetic Gravity Braiding (including the Galilean fluid) are conservative — they do not 
have a notion of entropy — but they are nonetheless imperfect fluids. One of the peculiarities of these 
systems is that the local properties not only depend on the fluid variables, which in the cosmological 
context are Lp and m, but also on the expansion (Hubble parameter) and on the external matter 
energy density p and pressure p, see equation of motion (3.11). In particular, the condition for the 
absence of ghosts (3.12) explicitly depends on H , i.e. on the external energy density, while the sound 
speed not only depends on the Hubble parameter, i.e. on p, but also directly on the external pressure 
p, because of the k term in (3.13) and the structure of the equation of motion (3.11). Thus, contrary 
to D, the sound speed also depends on the external equation of state w ~ pj p. This dependence 
of local properties of fluids with Kinetic Gravity Braiding on external p and p reveals a similarity of 
these fluids with open systems. 

These properties are crucial for our analysis and the appearance of the superluminality. 

4 Robust subluminality without external matter 

In this section, we will work with the n field and action (2.1) without any external matter. First of 
all, we eliminate 7 by rescaling (2.6) so that the time variable which enters both tt and H 

i = 7i/V, (4.1) 

and consequently the rescaled variables m and h are defined by 

TO = TT 7I/2 , and h = H 7I/2 . (4.2) 

We will denote differentiation with respect to this new time variable with a prime, ( )' = d/dr. 

Following our paper [36] it is useful to describe dynamics in phase space (to, h) by solving the 
Friedmann equation (3.9) with respect to e^'^: 



2^ _ 3to4/3 + 12/ito3 - 6h^ 



e 



2m2 



(4.3) 



Thus the expression on the r.h.s. must be strictly positive. Wherever the above function is negative, 
it corresponds to a region of the phase space which is dynamically inaccessible. The advantage of the 
variables {m, h) is that one avoids in this way branches in the square root and related folding of the 
phase space present for (tt, m) description, see discussion in the subsection 5.2. 

In the variables (to, h) the evolution of the system is given by equations (3.11) and (3.10) which 
after elimination of tt and rescaling can be written as 

, [l2/i2-m2 (to2- l)/3 + /(,TO (4- 6to2 + /3)1 to3- 2/1^(3/1 + to) ^ ^ 

TT^ — 1TI ^^ 7; V^ ?; r^ , (4.4) 

2/i2 + to4 (2to2 + /3) ' ^ ' 

, [I2h^m^ + (8/ito3 - 4/i2 + m^l3) (1 + /?) + to"/?] to^ - 12h^ 



2 (2/i2 + to4 (2to2 + /?)) 

We have plotted this phase space for a number of values of /3 in Fig. 1. 

This system evolves in the region <i>(TO, h) restricted because of (4.3) to: 

m'^P + Ahm^ - 2/1^ > . (4.5) 

Thus the phase space <I>(to,, h) is located between two curves 



h± (to,) = m 



m ± \l TO,^ H 



(4.6) 



Further, the absence of ghosts reqinres the positiveness of D from (3.12) and provides the bound 

-D ^ 2h^ + m" {2m^ + /3) > , (4.7) 



9 



which is always satisfied for /3 = 1 + a > 1. On Fig. 1, we have marked the inaccessible regions in 
deep red with the boundaries defined by (4.6). There are no ghosty regions. 
The general formula (3.13) for the sound speed translates for Cg (m, h) to 

Cs = 1 - 7-^ m^ {m, h) , (4.8) 

where the polynomial cubic in h 

^ (to, h) = 8h^ + 2h^m (l - im^ + /?) + 4/(,to* {2m^ - l) + m^ {'im^ + ^rn^l^ + (/3 - l)/3) , (4.9) 

defines whether the sound speed is less or larger than the speed of light. Subluminality corresponds 
to m^ (to, h) > 0. For each m, the polynomial ^ (to,, h) can have maximum 3 real roots hs (to). The 
number of roots, when solved for h, is defined by the sign of the discriminant & (a, to), 

^f = 2144to1° + 16to® fl25 + 203a]+2TO^ [339 + a(1066 + 915a)] + (4.10) 

32TO,'' 

+ 4 m'' [38 + a(126 + a(194 + 115a))] + to^ [32 + a(116 + a(154 + a(119 + 47a)))] + 
+ a(l + a)(2 + a)^. 

This is the only formula for which the parameter a = (3 — 1 is more convenient. The polynomial 
on the r.h.s. has only even powers of to and is manifestly positive for all non-negative a. Thus the 
discriminant ^ (a, to,) of the polynomial ^ {m, h) is always negative for all to, 7^ and non-negative 
a, whereas & (a, 0) — 0. Therefore, for to y^ 0, the polynomial 3^ (to, h) has only one real root hg (to,) 
which one can find analytically, but the expression for which is by far too long to be suitable for any 
analysis ^. Now we can evaluate the value of m^ (to, h± (to)) at the borders h± (m) of the phase 
space ^{m,h) given by (4.6). For ft,_ (to) we obtain 

TO^ (to, h^ (to)) = ^^ {e - ^e + 9^' - se + 4) , (4.11) 

(2to^ -I- p) 



where ^ — m/yjm? -1-/3/2. This function m^ {m, h- (to)) is always non-negative since the poly- 
nomial in ^ is positive for ^ = and has four complex roots ^ (2 — z ± \/— 5 — 4i) and the complex 
conjugates. For the other boundary, h+ (m), we have 

TO^ (to, h+ (to)) = ^J" {e + 4^' + 9^' + 8e + 4) , (4.12) 

(2to,^ + p) 

which is also always non- negative, because the polynomial in ^ has four complex roots of the form 
i (—2 — i ± ^/—5 — 4i) and the complex conjugates. 

This means that the single real root hs (to,) of £P (to, h) is not located in the phase space and 
m^ (to, h) > not only on the boundaries h± (to) but the whole cosmological phase space ^{m, h). 
Therefore, for all a > 0, we conclude that the system is suhluminal on all cosmological configurations. 
In this sense we have proven a strong version of the desired property 6 on page 2 from [54] . 

We have illustrated this in Fig. 1. The roots hs{m) of (4.9) are the boundaries of the blue regions. 
As can be seen for /3 > 1 (a > 0) there is only one such boundary and it always occurs inside the deep 
red dynamically inaccessible region, as we have proven. Only in the case of /3 = (a = —1), which 
is not part of the parameter space of Suhluminal Galilean Genesis, a new region of superluminality 
appears: now, for a range of values of to, there are three real roots hs{m) delineating the boundaries 
of these superluminal regions. Two of these roots lie inside the phase space. 



"Note that for a < this situation changes so that for sufficiently small m the polynomial .^ (m, h) has 3 real roots 
and there is always a superluminal region, see the lower right panel in Fig. 1. 
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5 Superluminality with external matter 

In this section we consider cosmology with the Galileon in the presence of external matter with a 
standard equation of state w = p/ p = const, e.g. dust or radiation. We will restrict our attention to 
matter satisfying the Strong Energy Condition (SEC) and the Dominant Energy Condition (DEC) 
so that —1/3 < w < 1. The main purpose of this section is to prove that when such external mater 
is added into the game there are some regions in cosmological phase space where the fluctuations of 
the Galileon propagate faster than light. Moreover, in some of these regions the perturbations are 
not ghosty so that these configurations should belong to the same effective filed theory (EFT) which 
describes the regions where the perturbations are subluminal and those where the external matter is 
negligible or absent at all. 

5.1 {m,h,p) coordinates and a proof of general superluminality 

First, similarly to the previous section, we can proceed by working in the phase space (m, h, p) where 
p is rescaled as in (2.7). One can eliminate vr by solving the first Friedmann equation (3.9) with 
respect to tt so that 



2p + 3m^P+l2hm^ -6h^ 



2m2 
Therefore the phase space is restricted by 

-p-2h^ +4:hm^ + 771^^13 >0. (5.1) 

o 

The perturbations of the Galileon are not ghosts when D from (3.12) is positive which provides the 
inequality 

m^D = 6h^ + 6m^ + 3m^/3 - 2p > . (5.2) 

We would like to stress that, contrary to k-essence, the coefficient D explicitly depends on the external 
energy density p. Thus a sufficiently large positive amount of the external energy density can turn 
the perturbations of the Galileon into ghosts. Moreover, in the presence of p there are always ghosts 
close enough to the origin of the phase space, m = h = 0. 

From the above inequalities, it follows that the external energy density should be located 

3h^ + -m^/3 ~ 6hm^ < p < 3h^ + -m^fi + 3m^ . (5.3) 

This is only possible for such (to, h) that 

TO (m^ + 2h) > . (5.4) 

In particular this works for to, > and /i > 0. For the sound speed (3.13) one obtains 

cl = l 2 ^ [^^ ("^' ^)- P (8^ + 2to (/? + 1) + TO^ (1 + 3w))] , (5.5) 

{Dm) 

where £^ (to,, h) is given by (4.9) and is the same polynomial as in the previous section. For m > Q 
and 8ft, + 2to, (/3 + 1) + w? (1 + 3^) > 0, the sound speed is larger than the speed of light provided 

3^(m,ft) 

P > Psup - 8/^ + 2m (1 + /3) + m^ (1 + 3w) ' ^ ' 

If, contrary to the matter-free case, dynamics (5.1) allow for such {ra^h) that m3^ {m,h) < 0, then 
a positive amount of external matter which satisfies the no-ghost condition (5.2) will make the per- 
turbations superluminal. 
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Now let us prove that it is always possible to find such p. It is convenient to introduce a "probe" 
energy density 

/9A(TO,/^) = 3/^2 + ^m'*/3 + Am^ (5.7) 

We restrict the numerical parameter A to < A < 3. This restriction ensures that probe energy 
density is always positive and fits into the borders (5.3) for tti > and h > Q. On this probe energy 
density we have 



2 _, , Z{m,h) 



^^'-^+W^lAi#' (^-^^ 



where 



Z (m, h) = 18/1^(3 + w) + 8/iTO [3(1 + /?) + 2m^{\ - 3)] + (5.9) 

+ m^ [1213 + m2(9(u; - 3)/3 + 4(1 + /?)A) + 2m^{X (1 + 3w) - 12)] . 

A = 3 is the boundary of the ghosty region on which the sound speed (5.5) diverges. For normal mat- 
ter with w > —3, the polynomial Z as a function of ft. is a convex parabola. Thus for sufficiently large 
h this function Z is always positive. On the other hand, for given (m, h), we can choose A arbitrarily 
close to 3, i.e. arbitrarily close to the ghosty region. Thus we conclude that the sound speed of the 
Galileon can not only be superluminal but can, indeed, acquire arbitrarily large positive values. We 
would like to stress that by construction the perturbations are not ghosts on these configurations. 
This conclusion only uses the positivity of /3 which is needed for the positivity of p\ . 

One may think that the superluminality can only correspond to large h but in fact there are 
always regions around h = where the sound speed is larger than 1. Indeed, in that case 

-m^P < p<Pgh= -m'^l3 + 3m^ , (5.10) 



and 



2 _ _^ 4m4 [3m* {Am* + 5m^l3 + (/3 - l)/3) - p (2 (/3 + 1) + m^ (1 + 3w))] 

(6m6 + 3to4/3 - 2p)^ 



Thus, for w > —1/3 which are considering here, one has to require that 

_ 3m* {Am* + hm^j3 + (/3 - l)/3) 
^ ^ ^^"P " 2{P+l) + m^l + 3w) ' ^^-^^^ 

and there are p realising this regime provided psup < Pgh- This happens for those m for which 

Qm*{w - 1) + m^ (4 - /3(5 - 3w)) + 4/3 > . (5.13) 

The latter inequality is manifestly satisfied for sufficiently small m? because /3 > 1. Thus the super- 
luminality without ghosts can also happen for small m and h or in other words for p <C 1 . 

5.2 {(f), m,p) coordinates and induced superluminality of original configurations 

The disadvantage of choosing the phase space in terms of {m,h,p) in section 5.1 is that adding 
external energy density p while keeping both m and h fixed forces the value of the scalar field tt to 
change. Thus, in some sense, the new configuration including p does not correspond to the initial one. 
To avoid this, in this section, we will consider the phase space in terms of the coordinates (0, m, p), 
where (j) is given by (2.3), m — tt' = cj)' /(f>, i.e. is still given by (4.2) and p is again rescaled as in (2.7). 
However, we will see this is much more complicated as a result of the fact that this three-dimensional 
phase space is double folded, because of the square root 

np = ^/6 [2p + to2 (3^m2 -I- 6m* - 2(/)2)] , (5.14) 
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in the solution of the Friedmann equation (3.9) 



h = m^ + -np. (5.15) 



Here f2p = f2 (0, m, p) is defined in such a way that the sign of flp changes when the argument of the 
square root evolves trough zero. Only regions where the argument of the square root in f2p is positive 
correspond to physically available configurations, which result in a condition on p, 

p > (j)^ - '^m^ {P + 2m^) . (5.16) 

In this discussion, we are only considering normal matter and parameter ranges relevant to the 
suhluminal Galilean Genesis model, i.e. we will assume in what follows that 

— 1/3 < w < 1 , w — const (5.17) 

1 </3<4. 

Let us introduce the configurations of the Galilean which are dynamically allowed without any 
external matter: 

3 

r = {(j)Trn) : such configurations that 0^ < -m^ {(3 + 2to^) . (5.18) 

These configurations are exactly those which we have studied in the section 4, just in different variables. 
In the rest of this section we will only consider these configurations F. If we add any positive p > 
to any of these configurations F, it is still dynamically allowed, since (5.16) holds. The addition 
of the positive energy density p opens up the phase space — allows one to probe new cosmological 
configurations {4>,m) violating the condition (5.18). 
For the positivity of p we have to require that 

p positive: fll > fll ^ 6m^ [3/3™^ + Qm* - 20^] . (5.19) 

On configurations F we always have fl^ > 0. 

As we have shown before, in the section 4 all configurations F have subluminal sound speed when 
there is no external matter. In the subsection 5.1 we have showed that superluminality takes place 
somewhere on available and not ghosty phase space. However, these regions could not be available 
without external matter - so that for {4>,m) the condition (5.18) is violated. Our main goal now is 
to prove that around some of the configurations F the sound speed can become superluminal on the 
addition of external matter, but the perturbations remain not ghosty. 

The perturbations are not ghosts when D from (3.12) is positive, which provides in our current 
variables the inequality 

<j)^D , , (/)2 - 3m2 (/3 + 3m2) 

— — = m (fip - righosts) > , where fighosts = ■ (5.20) 

Therefore the condition for the absence of ghosts reduces to 

No ghosts: fip > righosts , for m > , 

^p < f^ghosts , for m < , (5.21) 

Note that fighosts (r) > for m < and fighosts (r) < for to > 0. Further it is convenient to express 

r^^ + 18m4 (/3 + 4to2) 

"ghosts ^ Z~Z ^ . (0.2Z) 

iZm,'^ 
Hence, on F and for sufficiently small m^ one always obtains ^i^jj^g^g (F) > fip (F). 



13 



The sound speed (3.13) as a function on phase space {4>^m^p) is given by 

2, . ■m^n]{(t>,m,np) 



where 



^^ {(j), m, np) = (7 - 3w)nlm + WUp (3m2 (/3 + 3™^) - c/)^) + 
+ 6m [3m^ (4/3^ + 3(9 + w)Pm'^ + 6(7 + w)m'^) - 2(j)'^ (2(1 + 13) + 3(3 + w) m^)] 



(5.23) 



(5.24) 



The superluminaHty is present in regions where m^^^ {<f>, m, ilp) < 0. Mw (</>, rn, flp) depends on p 
only through flp, with respect to which m^w (0, ^nn, fi^) is a convex parabola. Thus superluniinality is 
only possible if the roots ^w± of !^.w (0, "^, ^p) are real. Given that, in order to have superluminality, 



ilp must lie between these real roots: 



Superluminality exists: 






< VI p < Vlw+ , 



(5.25) 



7-3w 



n 



ghosts 



=•= -^'^■5hosts + ^ 



, ^w± G 



where fighosts is defined in (5.20) and 



(7 - 3w) 
16 



9m2 [(/3 - l)/3 + (7/3 - 2)m2 + 12771*] + [2(1 + /?) + 3(3 + 70)777^ 



"0 

4777^ 



(5.26) 



with Q.Q given by (5.19). Now we will prove that there are such {(j),m) from F that these roots exist 
and that at least some fip between these roots do not have ghosts and have p > 0. 
For riu,± of (5.25) to be real, we need to require that 



f^ghosts+S>0. 



(5.27) 



For our chosen parameters w and (3, (5.17), S (F) < 0, on all configurations F. Therefore, if 51u,± 
exist somewhere on F, then for to, > we have i^w- (r) < ^w+ (r) < whereas for tti < these roots 
are located < ^10- (r) < flw+ (r). 

To avoid ghosts as in (5.21), we have to require that 



ii-U}+ ~ 



i ii/t — 



7-3it; 


f^ghosts + y ^ghosts + ^ 


8 
7-3i7; 


fighosts - y ^ghosts + ^ 



> f^ghosts , for m, > , 
< f^ghosts , for 777 < . 



(5.28) 



For configurations F and normal matter these are equivalent to the condition 

2 

j ( i + 611) ) 

ghosts 



Superluminal & no ghosts: T, > fl^ 



1 



(5.29) 



which is a stronger condition than the condition (5.27) for the existence of roots ^w± and which 
implies 

Si,„_ < s2grhosta < ^w+ ■ (5.30) 



Therefore if (5.29) holds and m^ is sufficiently small, one can always adjust ftp to be sufficiently 



close to righosts so that Qq < Ql < ^ly^^^t,, and therefore have superluminality with p > and no 



ghosts. 

Now let us find simple sufficient conditions on ((f), m) from F to satisfy (5.29), proving that there 
are accessible superluminal configurations. 
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The inequality (5.29) is equivalent to 

J"(02j = (3 + w)</>'' + S(/)2 + C>0, (5.31) 

where 

B = 2m^ [4 - (5 + 3w) /3 - 3m2(3 + w)\ , (5.32) 

and 

C = Zm" \im^{hw - 1) + Ylm^w^ + (1 + 3u;)/32] . (5.33) 

If for some range of m there are no real roots for the equation J^ (0^) — 0, then the inequality (5.31) 
(or equivalently (5.29)) holds for all </> and in particular for such (/> that the matter-free configuration 
belongs to F, (5.18). Inequality (5.31) holds for all 0, provided the quadratic equation for 0^ does 
not have positive solutions, in particular, if the discriminant is negative, i.e. if 

S{w?) = 18m^(l - w){i + w) - 3m^(3 + w)(4 + (3u; - 5) /3) - 4(5 + 3w - 2/3)/3 + 8 < . (5.34) 

Let us consider sufficiently small w? . For any /? > 1 and 

w>Wfi^ -^ '^ , (5.35) 

we obtain S (0) < 0, so that for sufficiently small vr? we have S (m^) < 0. In particular for /3 = 1 we 
have w\ = — i/s thus any SEC satisfying matter creates a small region close to to = such that the 
sound speed is superluminal and yet the perturbations are not ghosts. For /3 = 2 this translates into 
W2 — Q '&0 that reasonable matter would always create a superluminal region. Whereas, for /? = 4, 
which is the maximal interesting value of this parameter from [54] we have W4 = ^/e which would 
violate the DEC and which is not an equation of state available for standard matter. The limiting 
equation of state Wfi corresponds to radiation when 

Aad = ^^^ - 2.62 , (5.36) 

and to ultra-stiff equation of state Wfi = 1 when 

/3stiff' = 2 + 73 ~ 3.73 . (5.37) 

Thus for 1 < /? < /3stiff the addition of standard matter creates a superluminal but not ghosty region 
at least for those configurations from F which have sufficiently small m. 

Now we can look for more general larger to. Since S(vr?^ of (5.31) is a convex parabola, if we 
require that S{m'^) = has at least one positive root m^, then there will be a range of to where the 
inequality (5.34) is satisfied and where consequently (5.29) holds for all (j>. These roots are given by 



2 12 - 15/3 + wjA + (4 + 3w)(3) ± ^(3 + w)Y {w, /?) 

'm+ = TTTT^ TTT^ ^ , (5.38) 

^ 12{1 ~ w){3 + w) ' ^ ' 

where 

y (w, /3) = (1 + w)(ll + 3w{3w - 4))/?^ - 8(1 -I- w){9w - 5)/3 + 80w - 16 , (5.39) 

which is a quadratic polynomial in /3 and cubic in w. We now must ensure that both y > and 
rn^ > 0. For normal matter (1 -|- z<;)(ll -|- 3w{3w — 4)) > 0, thus for Y to be positive, the parameter 
/3 must be larger than 

4 U(4 -I- 9w) + 2^3(3 + w)(l-w2)(l-3w) - 5 
^* (") - (1 + ^(11 + 3^3.. -4)) ' ^'-''^ 
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for normal matter with w < 1/3, and /3 can be arbitrary for w > i/s- It is easy to check that /3* < 1 
for all w we consider, thus Y is always positive for our set of parameters (5.17). 

For < w < 1/3 one obtains that m'^ is always negative whereas m'^ is positive for either 



2 < /3 < Aad and ^^ ^^^f ^^ < w < ^- . (5.42) 



1 < /3 < 2 and < w < 1/3 , (5.41) 

or for 

3-1) 1 

3/3 

If we allow for matter with an ultra-relativistic equation of state 1/3 < w < 1 then 

Aad < /3 < /3stiff and ^^ ~ ^3 ^ " ^^ <w<\. (5.43) 

Thus we have proved that for any model with /3 < /3stiff — 3.73 one can add a positive amount of ex- 
ternal energy energy density with < w < 1 to make the sound speed of the perturbations around at 
least some original configurations larger than the speed of light. We presume that this statement also 
holds for /3stiff < /3 < 4. But to prove this it would require a more cumbersome analysis e.g. requiring 
the roots <l? of T (i/)^) = to be negative. We should also stress that we have only constructed one 
type of configurations with superluminality. We have not performed an exhaustive search to classify 
all such configurations. 
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